Spin-orbit effects in Na4lr308, a hyper-kagome lattice antiferromagnet 
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We consider spin-orbit coupling effects in Na4lr308, a material in which Ir'*''' spins form an 
hyper-kagome lattice, a three-dimensional network of corner-sharing triangles. We argue that both 
low temperature thermodynamic measurements and the impurity susceptibility induced by dilute 
substitution of Ti for Ir are suggestive of significant spin-orbit effects. Because of uncertainties 
in the crystal-field parameters, we consider two limits in which the spin-orbit coupling is either 
weak or strong compared to the non-cubic atomic splittings. A semi-microscopic calculation of 
the exchange Hamiltonian confirms that indeed large antisymmetric Dzyaloshinskii-Moriya (DM) 
and/or symmetric exchange anisotropy may be present. In the strong spin-orbit limit, the Ir-O-Ir 
superexchange contribution consists of unfrustrated strong symmetric exchange anisotropy, and we 
suggest that spin-liquid behavior is unlikely. In the weak spin-orbit limit, and for strong spin-orbit 
and direct Ir-Ir exchange, the Hamiltonian consists of Heisenberg and DM interactions. The DM 
coupling is parametrized by a three component DM vector (which must be determined empirically) . 
For a range of orientation of this vector, frustration is relieved and an ordered state occurs. For 
other orientations, even the classical ground states are very complex. We perform spin-wave and 
exact diagonalization calculations which suggest the persistence of a quantum spin liquid in the 
latter regime. Applications to Na4lr308 and broader implications are discussed. 

PACS numbers: 75.10.-b,75.10. Jm,75.25.+z 



I. INTRODUCTION 



Geometrically frustrated antiferromagnetism is a rich 
subject enjoying considerable theoretical and experimen- 
tal attention over several decades of researchi^^ Such 
systems are realized by materials containing magnetic 
ions in which the strongest antiferromagnetic exchanges 
occur on a network of bonds containing many trian- 
gular units. The most celebrated examples are the 
two-dimensional kagome (corner-sharing triangles) lat- 
tice and three-dimensional pyrochlore (corner-sharing 
tetrahedron) lattice. In ideal classical models, these lat- 
tices support highly degenerate ground states which pre- 
vent order down to very low temperature. Instead, the 
spins continue to fluctuate strongly despite significant 
correlations induced by the frustrated interactions. Sys- 
tems in this regime are dubbed (classical) spin liquids, 
or cooperative paramagnets. A major goal in the field 
is to ascertain whether such spin liquids might also oc- 
cur even in the zero temperature limit, in which both 
quantum effects and many non-ideal features of the ma- 
terials must be taken into account. The answer to this 
question is quite subtle, due to many competing effects 
that can come into play. Quantum and thermal fluctua- 
tions may break the ground state degeneracy and actu- 
ally induce magnetic order, an effect known as order-by- 
disorder i^i^i^ This effect, however, is understood theo- 
retically only in the large spin, S* 3> 1 limit, in which 
spins behave semi-classically. Nevertheless, some mod- 
els even with the smallest possible spins, S — 1/2, seem 
at least qualitatively to follow the order-by-disorder sce- 



nario. Conversely, in other models with small spin, quan- 
tum spin liquids have been shown to occur. No general 
theory to predict which of these two tendencies obtains 
exists at present. 

Despite this lack of theoretical discrimination, exper- 
imentalists have forged onward in recent years, uncov- 
ering a number of promising candidate quantum spin 
liquid materials with small spin S = 1/2 on geometri- 
cally frustrated lattices. These include an organic mag- 
net, K-(ET)2Cu2(CN)3, containing spin-1/2 moments 
on a slightly spatially anisotropic triangular lattice, 
ZnCu3(OII)gCl2, an inorganic realization of a spatially 
isotropic spin-1/2 kagome antiferromagnet, and very re- 
cently the cubic material Na4lr308 which realizes an 
hyper-kagome antiferromagnet, in which spin-1/2 mo- 
ments reside on a three-dimensional network of corner- 
sharing triangles^ - see Fig. [T] None of these compounds 
exhibit indications of magnetic ordering. The interpre- 
tation of the first two materials, however, is compli- 
cated by the appearance of inhomogeneous magnetic mo- 
ments at low temperature in k-(ET)2Cu2{C1SI)^, and by 
fairly high levels of substitutional disorder (Zn for Cu) in 
ZnCu3(OII)6Cl2. By contrast, the Ir''^ moments are ex- 
pected to be well ordered in Na4lr308 , due to the much 
larger ionic radius of Ir compared to Na and O. 

Recent two works^'^" assumed the nearest neighbor 
antiferromagnetic Heisenberg model for Na4lr308. In 
Ref. 9, the authors treated the spin as a classical 0{N) 
spin. By a large- iV mean field theory and classical Monte 
Carlo simulation, they found that the classical ground 
states are highly degenerate and a nematic order emerges 
at low temperatures in the Heisenberg model {N = 3) 
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via "order by disorder", representing the dominance of 
coplanar spin configuration. In Ref. [Tol . the authors 
presented a large-A^ Sp{N) method and studied both 
the semi-classical spin and quantum spin regimes. In 
the semi-classical limit, they predicted that an unusual 
k = (0,0,0) coplanar magnetically ordered ground state 
is stabilized with no local "weather vane" modes. While 
in the quantum limit, a gapped topological Z2 spin liquid 
emerges. 

Due to the large atomic number {Z = 77) of Ir, how- 
ever, we should carefully consider the role of spin-orbit 
coupling, whose leading effect in localized 5=1/2 elec- 
tron systems is the Dzyaloshinskii-Moriya (DM) interac- 
tion in the weak spin-orbit coupling limt. In fact, DM 
interactions have been argued to play an important role 
even in ZnCu3(OH)6Cl2, with much less relativistic Cu 
{Z = 29) momentsiii The DM interaction reduces the 
full SU(2) spin-rotational invariance of the Heisenberg 
Hamiltonian to the Z2 discrete time-reversal symmetry 
(in addition to coupling spin transformations to the dis- 
crete point group operations of the lattice). On gen- 
eral grounds, this is expected to lower the degeneracy of 
the classical ground state manifold. However, depend- 
ing upon the detailed form of the DM coupling, vary- 
ing degrees of degeneracy remain, indicative of different 
amounts of frustration. The tendency of the system to 
retain the order of the classical ground state is certainly 
also variable, and warrants investigation. This is one of 
the motivations of the present study. 

Another motivation comes directly from the experi- 
ments in RefJ., several aspects of which are suggestive 
of the presence of spin-orbit coupling. First, the "Wil- 
son ratio" R = Tx/cy, is observed to grow with coohng 
at low temperature, following a power-law R ~ 1/T"~^, 
with 2 < a < 3. Here x ^ const, is the magnetic suscep- 
tibility and Cy ^ T" is the specific heat. As will be dis- 
cussed in Sec. HIl such a low temperature behavior is in- 
compatible with any spin-rotationally invariant phase of 
matter supporting well-defined quasi-particle excitations. 
To our knowledge, it is at odds with all known theoret- 
ical models of quantum spin liquids, and seems highly 
unlikely on general grounds. Taking into account the ob- 
served field-independence (up to 12 Tesla) of the specific 
heat Cy brings the behavior even further into disagree- 
ment with spin-rotationally invariant theories. Second, 
samples in which a fraction x of Ir atoms are substituted 
by Ti (which are in a non-magnetic Ti^"*" state) display 
a Curie component in the susceptibility linearly propor- 
tional to X with a strongly suppressed amplitude, of ap- 
proximately 1/3 of a spin-1/2 moment per Ti. As we also 
show in Sec. Ill Bl such behavior is also at odds with any 
simple spin-rotationally invariant low-temperature phase 
(assuming no clustering of the Ti atoms), though some 
more exotic. All these observations, however, are read- 
ily reconciled by assuming the presence of spin-rotational 
symmetry breaking. Given the lack of any observed mag- 
netic ordering, explicit and substantial spin-orbit inter- 
actions would appear to be a likely candidate. 



In Sec. mil we consider an explicit semi-microscopic 
calculation of the exchange Hamiltonian in the presence 
of spin-orbit coupling. We consider both super-exchange 
through the intermediate O atoms, and direct exchange 
between closest pairs of Ir spins. The results depend cru- 
cially upon the relative magnitude of the spin orbit cou- 
pling constant A and the non-cubic splittings of the t2g 
multiplet. This is quantified by two dimensionless ratios 
of A to the two energy splittings £2 — ei and £3 — ei of the 
orbital levels in the absence of spin-orbit. When A is the 
largest energy scale - the "strong spin orbit limit" - the 
"spin" has a substantial orbital angular momentum com- 
ponent, while in the opposite "weak spin orbit limit", it 
is predominantly microscopic spin angular momentum. 
Indeed, the g-factor has opposite sign in the two lim- 
its. Which if either limit applies is the most fundamental 
physical question to be understood concerning the nature 
of magnetism in Na4lr308. We are not aware of any cal- 
culations or direct experimental measurements that indi- 
cate whether Na4lr308 is in the weak or strong spin-orbit 
limits, or intermediate between these situations. Instead 
we will address this question by comparing the expected 
phenomenology for the two cases to experimental obser- 
vations. 

In the strong spin-orbit limit, when the dominant 
mechanism is Ir-O-Ir superexchange, we find an highly 
anisotropic effective spin Hamiltonian, in which two spin 
components on each bond interact antiferromagnetically 
while the third interacts ferromagnetically. Specifically, 

H = ^J45f5j', (1) 

iv) 

where (ef^ , e^^ , ef^ ) is a permutation of (-1-1, +1, —1) cho- 
sen appropriately for each bond (see Sec. lIIIF]) to specify 
the two antiferromagnetic and one ferromagnetic direc- 
tion. We call Eq. ([T|) the "strong anisotropy" Hamilto- 
nian. 

Somewhat surprisingly, the remaining three cases: 
strong spin orbit and direct exchange, and weak spin- 
orbit and superexchange or direct exchange, all lead to 
approximately isotropic Heisenberg interactions. For the 
weak spin-orbit limit, this is guaranteed, but it is cer- 
tainly not in the strong spin-orbit case. The dominant 
spin-rotational symmetry breaking effect, which is per- 
turbative in all three regimes, is the DM interaction. The 
effective Hamiltonian has the form 

7i = ^[JS..S,+Dy.(S. xS,)]. (2) 

{ij) 

Here J is the same for all bonds, and estimated as 
J « 400ii' from the measured Curie- Weiss temperature 
©cw ~ —650K. Symmetry strongly restricts the struc- 
ture of this effective magnetic Hamiltonian for hyper- 
kagome. The full set of DM vectors D^^- may be fixed by 
just three parameters. That is, D^^ on any one bond is 
arbitrary (by symmetry), but that choice determines all 
remaining Dy in the system. It is convenient to choose 



FIG. 1: The hyper-kagome lattice of Ir*+ spins, with one 
classical ground state of the strong anisotropy Hamiltonian 
shown. Though this particular ground state is collinear, other 
ground states are not. 

the local coordinate system {Di, D2, D3), where Di is 
the component aligned with the bond, D2 is normal to 
the triangle plane in hyper-kagome lattice, D3 lies in the 
triangle plane but perpendicular to the bond (see Fig. [3]). 
The semi-microscopic calculations in Sec. IIII Gl confirm 
that all three components are non- vanishing, and gives 
a quantitative understanding of them. Due to the large 
A and considerable uncertainties in estimating the non- 
cubic energy splittings, it is difficult to estimate the over- 
all magnitude of the DM terms, but there is no reason 
they need be particularly small, though the pcrturbative 
estimates are presumably valid only for \Di\ < J < 0.1 
or so. A naive estimation is obtained by noting that 
in this limit the ratios of DM to exchange are expected 
to be of the same order as the shift of the g-factor, 
i.e. \Di\/J ~ |1 — |g|/2|. From the measured moment 
fj-eff ~ 1.9/i_B = IfflM-B, assuming we are in this limit 
would give \Di\/J ^ 0.05 or so. In Sec. IIVI wc consid- 
ered the strong anisotropy Hamiltonian, Eq. ^ in the 
classical approximation. Remarkably, unlike the Heisen- 
berg model which is macroscopically degenerate (i.e. its 
ground states are specified by a number of continuous 
parameters proportional to the number of spins), the 
system in this limit has an almost unique ground state. 
We find a continuous two parameter manifold of ground 
states, in which any one spin can be specified arbitrar- 
ily after which all others are determined. This is still a 
(small) accidental degeneracy, since the system has itself 
only discrete (space-group and time-reversal) symmetries 
which do not protect any continuous degeneracies. Nev- 
ertheless, this degeneracy is presumably insufhcient to 
prevent ordering in a classical system. The behavior in 
the physical S = 1/2 quantum problem is not known, 
but one would expect that an ordered phase of the same 
symmetry as the classical one is rather likely, and there is 
little reason to suppose a significant suppression of the or- 
dering temperature relative to the Curie- Weiss scale. The 
disagreement of these expectations with the experimen- 



FIG. 2; (Color online) The "windmill" state, which is the 
classical ground state in the weak anisotropy limit when D2 < 
0. It is also the basis vector ■i/)2 (Table. [V]! of one dimensional 
representation F^'' (see Eq. (|62|l '). In the generic system with 
non-zero _Di , , the spins are slightly canted out of the plane 
of each triangle. 



tal observations suggests that it is the weakly-anisotropic 
DM Hamiltonian rather than this one which is most ap- 
propriate. We however return to this question in more 
detail in Sec. IVlIl 

In Sec. |Vl we turn to the weak anisotropy limit, and 
first explore the classical phase diagram of Eq. In 
general, even this optimization problem is highly non- 
trivial, given the large unit cell of the hyper-kagome lat- 
tice, and the possibility that the magnetic unit cell of 
the ground states may be yet larger. In the special case 
Di = D3 = and D2 < 0, however, it is possible to 
solve this problem exactly. The degeneracy is broken 
completely to a single Kramer's pair of coplanar ground 
states, for which the magnetic unit cell is equal to the 
crystallographic one. These may in this sense be con- 
sidered k — (0, 0, 0) states. One is drawn in Fig[2l We 
call this the "windmill" state. By several approximate 
methods, we establish the form of the phase diagram in 
the general D1-D2-D3 parameter space. Generically the 
windmill state distorts to a "canted windmill" state (still 
with k ~ (0, 0, 0)), occupying a finite region of the phase 
diagram. In addition, one finds a wide range of incom- 
mensurate phase, in which the ordering wavevector k is 
non-zero and generically irrational in reciprocal lattice 
coordinates. Owing to the breaking of space group sym- 
metries, the incommensurate phase retains more of the 
frustration-induced degeneracy. 

A key question is whether the DM interactions, ex- 
pected on physical grounds and invoked phenomenolog- 
ically to explain the experimental properties discussed 
above, are consistent with the observed spin liquid behav- 
ior of Na4lr308, i.e. the lack of any ordering down to the 
very low temperatures of T « 1.8K = Ociv/360. The 
breaking of degeneracy by DM might be expected to re- 



4 



duce quantum fluctuations and thereby lead to ordering, 
in conflict with experiment. To study this possibility, we 
carried out spin wave calculations of the excitation gap 
and the quantum correction to the classical ordered mo- 
ment. Indeed, we find that deep inside the k = (0, 0, 0) 
phases, the quantum correction is not too large, which 
leads us to expect that the spin- 1/2 system exhibits the 
classical order. However, we find very large quantum cor- 
rections elsewhere in the phase diagram, even for fairly 
substantial \Di\. In our results, small excitation gap will 
lead to a large quantum correction to classical ordered 
moment. Decreasing the excitation gap by changing the 
DM vector will eventually destroy the classical ordered 
moment completely. In this regime, the large quantum ef- 
fects invalidate the spin- wave treatment and indeed leave 
open the possibility of a quantum spin liquid, consis- 
tent with experiment. To further confirm the results and 
treatment of spin wave theory, we implemented exact di- 
agonalization on a small cluster (six triangles with 13 
spins). The excitation gap obtained from numerical data 
of specific heat qualitatively agrees with the prediction 
of spin wave theory. 

The remainder of this paper is organized as follows. 
In Sec. mil we discuss the symmetry allowed DM vec- 
tor components and calculate the exchange spin Hamil- 
tonian with a microscopic theory for both strong and 
weak spin-orbit coupling. In Sec. IIVI we discuss the clas- 
sical ground states of the strong anisotropic exchange 
Hamiltonian obtained from Ir-O-Ir superexchange in the 
strong spin-orbit coupling limit. In Sec. |V] we will turn 
to look at the weak anisotropy Hamiltonian, namely, the 
nearest-neighbor Heisenberg model with small DM in- 
teractions. We first present the magnetic ordered state 
when I?2 < then discuss the more general case when 
nonvanishing Di and D3 components are present in the 
system. In Sec. I VII we present a linear spin wave theory 
to find the zero temperature quantum correction to the 
magnetically ordered phase and compare with exact di- 
agonalization. Finally, a discussion of our main results 
and their relevance to Na4lr308 is given in Sec. IVIII 



II. THERMODYNAMICS OF 
SPIN-ROTATIONALLY INVARIANT MAGNETIC 
PHASES 

In this section, we discuss some apparent constraints 
on the low temperature susceptibility and specific heat 
in spin-rotationally invariant phases of matter. As de- 
scribed in the introduction, these constraints appear to 
be violated in Na4lr308, which we take as an indication 
of the presence of substantial spin-orbit interactions. 



A. Clean system 

We take spin-rotational invariance to mean the exis- 
tence of global SU(2) spin symmetry. According to stan- 



dard quantum mechanics, this implies that all states may 
be chosen as eigenstates of ^^'^ "^tot' where S'tot 
is the operator for total spin. The choice of z axis being 
arbitrary, we take it along the axis of any applied field. 
The effect of the field on the system is then entirely de- 
scribed by the term 

Hh = -H'^SI = -HS^Qrj, , (3) 

i 

where we have absorbed the (presumed known) g-factor, 
Bohr magneton, etc. into the definition of H. One ob- 
serves from Eq. ^ that Hh is diagonal in the 
basis, and thus the Hamiltonian eigenstates themselves 
are independent of field, and only the eigenvalues change. 
Focusing on the states rather than their energies, we may 
say that the only effect of the field upon the system in 
equilibrium is to modify the occupation probabilities of 
states. In this sense, the magnetic field is a thermo- 
dynamic perturbation, and the susceptibility is a ther- 
modynamic quantity, determined only by the density of 
states. The specific heat is of course also such a thermo- 
dynamic quantity, determined from the same density of 
states. Thus they are connected. 

Specifically, the specific heat is a measure of the full 
density of states for all excitations above the ground 
state, irrespective of their spin quantum numbers. The 
susceptibility, however, only counts those excitations 
which carry non-zero spin along the field. The pos- 
sibility of spin-less excitations allows some independence 
of the two: by introducing more — states, one can 
increase arbitrarily while leaving x unchanged. How- 
ever, the converse is not true. It would seem difficult 
to increase x without also contributing to Cy. The only 
way in which this can be done is to introduce states with 
very large S'^ (which then contribute a large amount to x) 
but very low energy (and hence do not contribute much 
to Cy). This case corresponds to a system on the verge 
of a ferromagnetic instability. 

Without fine-tuning to such a point, we are led to ex- 
pect that, in the presence of SU(2) symmetry, the Wilson 
ratio, 

R=^ (4) 

should have an upper bound, corresponding to all excita- 
tions contributing both to x ^^'^ c^- This can indeed be 
shown provided we assume the system can be described 
by a non-magnetic ground state and non-interacting 
quasiparticles characterized by a spin 5'^ quantum num- 
ber. We define the density of state g^{e) and ^^(e) for 
boson or fermion excitations carrying spin — m, re- 
spectively. The specific heat is 

/•oo 

Cv^drJ^ dee[g^{e)nb{e) + gl,{e)nf{e)] , (5) 
where 

nb/f{e) = . (6) 
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One obtains 
klT 



sinh''(a;/2) cosh^(x/2) 



Now consider the susceptibility 

/■oo 

X = <9h^ / dem[g'^{e)nh{e - Hm) 

m JO 

One finds 



(7) 

(8) 



H=0 



X 



oo 

m" I dx 





9m 

[kBTx) gl^{kBTx) 



sinh^(a;/2) cosh^(x/2) 



(9) 



In the low temperature limit, we may approximate 
gU^ [ksTx) by its small argument behavior, which is usu- 
ally a power-law form: 

(10) 

One needs obviously > — 1 for the density of states 
to be integrable (and hence the cumulative distribution 
well defined). We will encounter problems with Eq. ^ 
if 7m ^ 1 foi' ^'^y w ^ 0. This could be fixed by the 
inclusion of a chemical potential, whose temperature de- 
pendence we have ignored, and as usual is necessary to 
avoid Bose condensation of free bosons at low T when 
their energy is close to zero. This effect, however, does 
not change any of the results, so we have excluded it for 
simplicity here. 

Given Eq. PH)) . the specific heat will be controlled at 



low T by the minimum exponent over all 7; 



b/f. 



7o — mm < 7, 



in 1^ 



One has 



(11) 



(12) 



with some constant Aq. The susceptibility is controlled 
by the minimum exponent for m 7^ 0: 



71 = min|7^/-''; m 7^ o| 
Note that by definition, 70 < 71. Then 

Then the Wilson ratio becomes 

R ^ RqT^ , 
and 



where Rq 



Aok' 



71 - 70 > . 



(13) 
(14) 
(15) 

(16) 



Because T > 0, the Wilson ratio cannot diverge on low- 
ering T, and unless T = 0, actually vanishes as T ^ 0. 
In defining the Wilson ratio, we have considered only 
the zero field specific heat. In a field, contributions from 
all excitations with m ^ will be field dependent. So 
unless the m = mode is dominant in c„, the specific 
heat should be expected to be field dependent. Con- 
versely, field independence of the specific heat requires 
that the m — excitations dominate Cv In this case, 
we have T > 0, and the equality is not satisfied. Thus 
a field-independent low-temperature specific heat would 
be expected to correspond to a vanishing Wilson ratio 
as T — > 0. This makes the observed divergence of R 
on lowering T in Na4lr308 even more at odds with the 
theoretical expectation for an SU(2) invariant system. 

A few comments are in order. First, while we have 
assumed power-law forms for the low energy density of 
states, this is not essential. We believe the lack of low 
temperature divergence in R(T) is very robust within 
the quasiparticle picture. Beyond the quasiparticle ap- 
proximation, the situation is less clear, and we do not 
have a definitive proof of this behavior of R{T). How- 
ever, we do not know of any single theoretical counter- 
examples in the literature for SU(2) invariant low tem- 
perature phases. 

If SU(2) symmetry (or more specifically, invariance un- 
der spin rotations about the measurement axis) is bro- 
ken, however, one readily and indeed almost generically 
observes this behavior. This is quite familiar from the 
case of ordered antiferromagnets in two or three dimen- 
sions. These are well-known to display a non-vanishing 
constant zero temperature uniform susceptibility xo and 
a power-law specific heat c„ ^ AT'^ due to spin wave 
excitations, hence a Wilson ratio obeying Eq. ([TS]) with 
however T = 1 — d < 0. This arises because the ground 
state itself is modified continuously by the introduction 
of a magnetic field. Semi-classically, the magnetic field 
leads to a smooth canting of the antiferromagnetic mo- 
ments in the field direction, linearly proportional to the 
applied field. 

This phenomena is, however, not limited to systems 
with spontaneous symmetry breaking. It occurs when- 
ever the effective Hamiltonian for the low temperature 
phase does not conserve the spin component along the 
magnetic field. As an extreme example, one may con- 
sider the case of two spin- 1/2 spins coupled together by 
antiferromagnetic exchange and DM interaction: 

- JSi • S2 - i^z • Si X S2 - i/ {S^ + S^) , (17) 

where we have chosen the DM vector along the z axis, 
and therefore oriented the field along x so that it couples 
to a non-conserved magnetization. One can readily diag- 
onalize the Hamiltonian, and find that in zero field has a 
unique ground state with a gap A = 1/2{J + ^,P + D'^). 
Nevertheless, the susceptibility is non-zero when D ^ Q: 



dSf 



X = 



dH 



H=0 



J (V J2 +D^ + J) 



(18) 
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Because of the gap, the specific heat of the dimer is acti- 
vated at low temperature, and hence the dimer's Wilson 
ratio diverges exponentially at low temperature. In gen- 
eral, a non-zero limit for the low temperature susceptibil- 
ity is always to be expected once SU(2) symmetry break- 
ing perturbations are taken into account. The specific 
heat, however, is insensitive to symmetry, and remains a 
true probe of low energy modes. 



B. Impurity susceptibility 

In Na4lr308, the introduction of non-magnetic impu- 
rities (substitution of Ti'*+ for Ir*+) was observed to give 
rise to a Curie component with a reduced effective mo- 
ment of fieff ~ (2/^5 )/3 per Ti. We would like to argue 
that a spin liquid state with such a large reduction from 
the moment of a free spin, 2//^, is unlikely in the absence 
of spin-orbit interactions, but quite likely when they are 
invoked. 

Suppose the Hamiltonian has global SU(2) spin- 
rotational symmetry in the absence of an applied mag- 
netic field. Then a spin-liquid ground state, which, by 
definition, does not break SU(2) symmetry, must be a 
spin singlet, i.e. a state of total spin S = 0. Its exci- 
tations can therefore by characterized by spin quantum 
numbers. Representations of SU(2) always have integer 
or half-integer spin, and in particular for all these the 
projection of the total spin along any field axis is a mul- 
tiple of 1/2. 

Now consider a single impurity. It may be a strong 
perturbation locally, but does not perturb the Hamilto- 
nian far from itself. Again presuming spin-orbit can be 
neglected, the ground state of this system should be a 
spin eigenstate, though not necessarily non-zero. Nev- 
ertheless, it can be classified by a total spin which is a 
multiple of an half integer. It is natural to expect that the 
ground state multiplet of a single impurity controls the 
impurity susceptibility (but see below). Allowing now for 
an external field, this is simply described as in the pre- 
vious subsection by Eq. Since the low energy states 
are still good representations of SU(2), and only the total 
spin projection enters Eq. ([3]), we will obtain an effective 
moment which is at a minimum (if it is non-zero) 2fiB 
per impurity. 

The caveat in this argument is the possibility of a 
Kondo-like effect. If the spin liquid state is gapless, 
then there is a possibility for an impurity moment to 
be "screened" by the bulk degrees of freedom. Still, the 
possibility of a fractional impurity moment is delicate. 
Most Kondo effects either completely screen the moment 
(as in the single channel case, leading to /ie// = 0) or 
to weaker temperature dependence of the impurity sus- 
ceptibility (e.g. Ximp ~ |lnT| in the two-channel model, 
which has a non-trivial Kondo fixed point). Thus most 
types of Kondo effect do not allow for such behavior. 
Recently, it has been suggested that some spin liquids 
might sustain a critical fixed line of Kondo fixed points. 



connected to the free impurity fixed point. This situa- 
tion can in fact lead to a renormalized Curie constant 
It would indeed be appealing should such an exotic pos- 
sibility be realized in Na4lr308, but we should allow for 
simpler explanations. 

As is well-known, the effective moment of ions in solids 
varies widely from the quantized values expected from 
SU(2) symmetric considerations. This is of course due to 
spin-orbit coupling. In general, with spin-orbit interac- 
tions present, the ground state of an impurity can be ex- 
pected to be a Kramer's singlet or a Kramer's doublet. In 
the latter case, it will behave energetically (i.e. in specific 
heat) as a spin-1/2 spin, but will have in general a non- 
trivial g-tensor describing its coupling to a field. This 
refiects a change in the effective moment. Thus there is 
no "quantization" of the effective moment once spin-orbit 
coupling is substantial. The observed fractional effective 
moment in Na4lr308 is perhaps another indication in this 
direction. 



III. SPIN-ORBIT COUPLING IN THE 
HYPER-KAGOME LATTICE 

In this section, we discuss the form of the spin-orbit 
modifications to the isotropic Heisenberg Hamiltonian. 
This is not directly calculable from semi-microscopic con- 
siderations without some assumptions about the local 
energetics due to crystal field splittings. Therefore we 
consider below a number of cases. 



A. Symmetry allowed DM vector components 

In several cases, we will find that the dominant ef- 
fect of spin-orbit coupling is to induce Dzyaloshinskii- 
Moriya (DM) interactions between the nearest-neighbor 
spins. Therefore before attempting any calculations, it 
is instructive to first consider the symmetry constraints 
upon them. Generally, DM interactions are rather highly 
constrained. For instance, they are absent if there is an 
inversion center between the two spins in question (this 
is not the case in Na4lr308). The compound Na4lr308 
has cubic symmetry, described by the space group P4i32, 
and consequently has a number of point group symme- 
tries. For our purposes, it is useful to consider a uncon- 
ventional set of generators of these symmetries. Specifi- 
cally, the full point group can be generated from the set 
of 180° rotations around a local C2 axis at each site. Due 
to this symmetry, all the hyper-kagome sites and bonds 
are equivalent. In Table. |Vl we list the directions of the 
C2 axes (-01 ) for every site in the unit cell (see Fig. [3] for 
the labeling). The C2 rotational symmetries relate the 
DM vectors of any two bonds. That is, given the DM 
vector on any one hyper-kagome bond, all others are de- 
termined. This one DM vector, however, is itself entirely 
unconstrained by the P4i32 symmetry. 

Since any single bond of the hyper-kagome is uniquely 
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associated with one triangle, it is natural to adopt a lo- 
cal coordinate system based on this triangle to describe 
the DM vector's components. We denote the compo- 
nent aligned with the bond Di, the component normal 
to the triangle plane D2 and the component normal to 
the bond but localized in the triangle plane D3. Three 
components have been illustrated in Fig. [31 If we select 
the direction of Di component axis by assigning a di- 
rection to one bond (arrows in Fig. [3]), the C2 rotation 
symmetry can generate the equivalent Di axis for other 
bonds (see Fig. [3]). In every triangle, there is a chirality 
of the Di axis of three edges, which can be considered as 
the direction of D2 axis. The cross product of Di and 
D2 directional vector generates the direction of D3 axis. 





FIG. 4: (Color online) Ir^^ and octahedron O"^ environment 
(Thin black line). Two neighboring Ir''"'" are denoted by A 
and B (In orange). A/B's six O^" are labelled as 1/1', 2/2', 
3/3', 4/4', 5/5' and 6/6' (In pink), in which, 2 and 3', 5 
and 6' label the same points. The distances between Ir*^ 
and O^" order this way: \A5\ = \A6\ = \B5'\ = \B6'\ > 
\A3\ = \A4\ = |B3'| = \B4'\ > \A1\ = \A2\ = \B1'\ = \B2'\. 
The C2 axis (thick dash line) orients along -^(1,-1,0) at 
Ir''+ A and ^(0, 1, 1) at Ir''+ B. Mapped to the ideal hyper- 
kagome lattice, A and B correspond to point 4 and 8 in Fig.|3] 
respectively. 



FIG. 3: (Color online) Left: One unit cell of the hyper-kagome 
lattice. The pink balls are occupied by magnetic ions, which 
are connected by dark black bonds. There are 12 sites in one 
primitive unit cell. The arrow from site i to site j corresponds 
to T)ij ■ (Si X Sj) in the Hamiltonian. We will call these arrows 
DM interaction path. Right: DM vector components illus- 
trated on one triangle. Di is the component which is aligned 
with the DM interaction path (Left). D2 is the component 
normal to triange plane. The direction is decided by the chi- 
rality of bond direction. D3 is the component perpendicular 
to the bond but in the triangle plane. 



state 


5d orbitals at A 


5d orbitals at B 


parity 


1 1) 


xy 


yz 


even 


|2) 


■j^{xz-yz) 


■^{yx + zx) 


odd 


|3) 


■^(xz + yz) 


-^{yx-zx) 


even 


|4) 


2 2 
X -y 


2 2 
y -z 


odd 


|5) 


2 2 
6z — r 


ix^ - 


even 



TABLE I: The parity sectors of 5d electron orbitals by C2 
rotation 



Such a parametrization may be applied not only for 
the hyper-kagome lattice, but for any lattice consisting 
of corner-sharing triangles, such as the slightly distorted 
kagome lattice of Fe/Cr-jarositesJ^'i^ii^ In that example, 
the Di component is forbidden by a mirror plane sym- 
metry. In Na4lr308, there are as we said no constraints 
on the Di, and we might naively expect all three com- 
ponents to be non-vanishing and comparable. We will 
investigate this by microscopic calculations below. 

B. Local electron energetics of Ir ion 

Before moving to the microscopic theory of spin-orbit 
interactions, we need to understand the electron energy 
levels of the Ir*+ ions. With coordinates taken from Ta- 
ble. I in Ref. 0, two Ir^+ and their surrounding are 



drawn in Fig. [H For A ion, the C2 axis orients along 
-^(1,-1,0). Under this symmetry operation, x — > —y, 
y — > —x and z — > —2. Accordingly, we can group the 
bd orbitals into even and odd parity sectors, as shown in 
Table. E 

A large cubic crystal field splits the Cg and t2g states. 
The surrounding octahedron is slightly distorted to 
further split all the three t2g states. Ultimately no de- 
generacy is protected because the C2 symmetry has only 
one dimensional irreducible representations. The ener- 
getic ordering of orbitals shown in Fig. [5] was determined 
by looking at Coulomb interaction from surrounding O^^ 
and ignoring the spin-orbit interaction. 
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5d orbitals of 
Ir ion 



t2g 



I3> 
I2> 



Crystal field 
splitting 



Lattice distortion 



FIG. 5: The splitting and electron occupation of 5d orbitals of 
ions in the absence of spin-orbit interaction. The states 
are defined in Table. U 



Microscopic theory of exchange spin 
Hamiltonian 



Though symmetry determines the allowed non-zero 
components of the DM interaction, it does not 
give any guidance as to their relative and absolute 
magnitudesji^iii^iii In this part, we will derive the ex- 
change spin Hamiltonian from a microscopic point of 
view and obtain expressions from which crude estimates 
of the magnitude of various terms can be obtaine d^^'^^i^^ 
We consider both the hopping between Ir and O orbitals, 
and direct hopping between Ir orbitals. We also assume 
that the eg-t2g splitting is much greater than the split- 
tings among the three t2g states so that we can com- 
pletely project out the two e^ states. The model is then 
of five electrons on on the t2g orbitals of every Ir^'^. Fol- 
lowing some notations in Ref.[l3, we can write the Hamil- 
tonian of the Ir and O sublattice as 



where, 



7i = Wo + + Ti-LS , 



(19) 



kn<T 



Ud 



9 '^jma"'jm'c7'"'i'm'a'0'i 



jmm'aa' 



Up 
2 



2 PknaPkn'cr'Phn'cr'Pkna, (20) 

knn' aa' 

jma k{j)n 
Uj') mm' 

Hls = A^£, -s, . (22) 

3 

k{j) denotes the of the neighboring site j, 

'^Ima creation operator of an electron with spin 



a of the mth 5d orbital of ith Ir ion, is the energy 
of this orbital, m will take 1,2,3. plna 
ation operator of an electron on the 2p„ orbital with 
spin a. The energies are measured from the lowest en- 
ergy level of the Ir 5d orbitals, and Ud and Up are the 
Coulomb interaction constants between holes on the Ir^"*" 
site and site, respectively. We assume that Ud and 
Up are orbital-independent and ingore other "Kanamori 
parameters" >^ the inter-orbital exchange coupling and 
the pair-hopping amplitude, which should be small com- 
pared with Coulomb interaction. We also ignore the 
Coulomb interaction between two eletrons on different in- 
termediate ions. Here tjm,kn denotes the transfer of 
an electron between the mth orbital of Ir^"*" ion j and one 
of the 2pn orbitals of the neighboring ions k. Simi- 
larly, tjj^j,^, is the matrix element for electron transfer 
between m and m' orbitals on two nearest-neighbor (in 
the hyper-kagome sense) Ir atoms. £j and Sj denote the 
orbital and spin angular momenta at the jih Ir'*+ ion, 
respectively, and A is the spin-orbit coupling constant of 
the Ir^'^ ion. 

In order to understand the electron occupation on 
each site, we collect the quadratic terms for each site 
in Eq. (|22p and write down the onsite Hamiltonian as 



mm' aa' 



with 



M 



(i) 



(23) 



(24) 



where cr is the Pauli matrix and liimm' is the matrix 
element of between the mth and m'th orbital of the 
ith Ir'^'^ ion. It is useful to note^^ that the vector of three- 
dimensional matrix orbital angular momentum operators 
projected into the t2g manifold is actually proportional to 
the vector of orbital angular momentum operators for the 
3 ordinary {px,Py,Pz) states, but with a proportionality 
constant of — l! That is, suppressing the m,m' indices. 



(25) 



where L is a canonical angular momentum operator with 
= -|- 1) = 2. This effectively makes the spin-orbit 
coupling term directly analogous to the familiar one from 
an isolated atom with spherical symmetry in ap shell, but 
with the sign of the spin-orbit coupling reversed. 



D. Strong and weak spin orbit limits 

Obviously the nature of the "spin" itself (i.e. the 
Kramer's doublet ground state of the single whole in this 
multiplet) is crucially dependent upon the strength of the 
spin-orbit interaction, relative to the non-cubic splittings 
£3 — £3 ~ ci- This determines the nature of the wave- 
functions of the Kramer's pair, for instance the degree 
to which the "spin" carries true electron spin angular 
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moment or instead orbital angular momentum. This is 
more fundamental than the exchange interaction, so we 
consider it first. 

1. Strong spin orbit 

In the strong spin-orbit limit, we can to a first ap- 
proximation ignore the non-cubic splittings, and we have 
simply 

=A£,-s, = -AL, -s,. (26) 

This is of course diagonalized by constructing eigenstates 
of the "total angular momentum" 

J, = L, + s,. (27) 

Because of the minus sign in Eq. the highest en- 

ergy doublet is simply the j = 1/2 Kramer's pair. This 
describes the wavefunction of the half-filled orbital. It is 
natural to define the effective spin operator in this case 
as 

S. =J.. (28) 

Clearly it is a strong mix of orbital and spin components. 
According to the Wigner-Eckart theorem, the matrix el- 
ements of Si, Li and Ji are all proportional. This enables 
one, with a little Clebsch-Gordan algebra, to arrive at 
an expression for the magnetic moment operator (in the 
j = 1/2 manifold) 

M, = -^lB{^^ + = +2/iBS,, (29) 

where /is is the Bohr magneton. Interestingly, this is the 
same magnitude but opposite sign as for a free electron! 
It will of course suffer corrections perturbative in (e^ — 
£3) /A, as one moves away from the strong spin orbit limit. 

2. Weak spin orbit 

Now consider the weak spin orbit limit. In this case, 
for A = 0, the half- filled doublet is simply the m — 2> 
orbital, with two possible "true" spin orientations. Thus 
we approximately have 

S, «s, + 0(A/(e, -e,)). (30) 

Now there is essentially no orbital angular momentum 
component to the spin (€i ~ 0), and one obtains 

M, = -2^iBS,{l + 0{- :)). (31) 

|ei,2 - esl 

Note the important sign difference from Eq. (129]) . This 
is the most fundamental physical distinction between the 
weak and strong spin-orbit limits. However, the magni- 
tude of the proportionality between the magnetization 
and spin - the ^-factor - is the same in both cases. 
This means that the simplest experimental measure, the 
Curie susceptibility, cannot distinguish the two possibil- 
ities. We will consider both cases below. 



E. General exchange formulation 

We now turn to the exchange calculations. Let us con- 
sider the general case first. We must deal with M'^^\ 
which is a 6 X 6 matrix. Diagonalize A^^*^ so that 
7Vf(') — Here, £ is a site-independent eigen- 

value matrix, and T^'^ is a unitary eigenvector matrix. 
A^(*^ has three different eigenvalues £1, £2 and £3, each 
has a two-fold degeneracy due to Kramers' degeneracy 
theorem. The effective spin operator Si will be defined 
to act in this doublet. In the strong and weak spin-orbit 
limits, we have exphcitly Eq. (^5]) and Eq. ([50]) . respec- 
tively. Furthermore, we define a new set of electron cre- 
ation and annihilation operators 

with a.ima annihilates an electron on the £,n state with 
spin a at site i. 

Without losing any generality, we assume that £3 > 
£1.2, then £1^2 states are fully occupied and £3 state is 
half-occupied, leading to a total spin-i at every site. 
Accordingly, the magnetic momentum operator (M^ — 
—IJ-B{f-i + 2si) at each site should be projected onto the 
Kramers' doublet ground states: 

Mb ^„ 

mnap 

= ^^iLspf^lSa ■ Si (33) 

with cr the vector of Pauli matrices. Also, GiJ js and the 
effective spin operator Si are defined as 

mnaf3 

S» = X! \43a'^apa^3(3 , (34) 

a, 13 

and Pi is the ground state projection operator: 

Pi = a^3^ I 0) I + a,3| I 0) ((/) I 0^3^^ . (35) 

Here | is the £1,2.3 fully-occupied state. In the last 

step Eq. ((33|) . J2cr '^ia'^i^o' ~ 1 been used. 

Let's go back to Eq. (f^^ . and express the micro- 
scopic Hamiltonian in terms of ajma and a^„g.- Given 
the Hamiltonian in Eq. (I19|) . which includes the largest 
Coulomb energy U but neglects the smaller Hunds-rule 
exchange coupling between electrons in different orbitals 
on the same atom (and other similar interactions), only 
hopping through the half-filled orbital contributes to the 
super-exchange interaction. This is in accord with the 
"Goodenough-Kanamori" rules, which state that the ex- 
change coupling contributed from a half-occupied orbital 
and a fully-occupied orbitals is much weaker than the 
one from two half-occupied orbitals. Thus, we only need 
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to focus on the hopping between the £3 orbitals, as half- as 
occupied orbital. The microscopic Hamiltonian is written 



-W — \ " t I \ " ft 1 \ " c t J. El \ " t t 

kncr knn'aa' jma jmm'aa' 

jk(j)n af3 {jj') a/3 



with 



and 



ma 



^j,kn — Y^ ')^jm,knT^2.m0'^0<^ ' 



2 

m,Q(/3 



''j3,j 3 / ^ 2 i^'i'^ 3a,m(7 m'cr,3a 

mm' ,acT 

mm' ,(T,aP 

where cr is vector of the three Pauli matrices. Now we may follow the standard perturbative treatment of superex- 
change. We consider separately the superexchange through the intermediate O^^ ions, and the direct exchange 
contributions. 



1. Superexchange through oxygen ions 

In this case the leading contribution is fourth order in hopping, i.e. a result of fourth order degenerate perturbation 
theory. We must include four "hops" between Ir^+ and ions, which consist of "hops" described by spin-isotropic t 
matrix elements, and "hops" given by anisotropic C matrix elements. One thereby obtains the exchange Hamiltonian 
as 

= Y^JSr ■ S, + Dy • (S, X S,) + S, ■ 7y • S,] (39) 

iv) 

with the first two terms the Heisenberg and DM interactions precisely as in Eq. ^ , and the third term the anisotropic 
exchange. The explicit formulae for the coupling constants are: 

J 4 ^ ^ ^ij.kn9kn,k'n' ^ji.k'n' (40) 

kn,k'n' 

"^ij 4^ ^ ^ ij .knQkn^k' n' ji^k' n' ^ij^knSkn.k'n'^ji.k'n'^ (41) 

kn.k'n' 

^ ij 4 ^ ^ [( V ij^knSkn^k'n' ^ ji,k'n' ji,kn9kn,k' n' ^ ij,k'n') f i^ij,kn ' 9kn,k'n' ^ji,k'n')] ■ (42) 

kn,k'n' 

The vector with arrow ^ or — > indicates that inner product is taken with the spin operator put in the direction of 
< — > 

the arrow. 1 is a 3 x 3 unit matrix. Sij,kn, ^ij,kn and gkn.k'n' are given by 

^ij,kn — ti3,knikn,j3 ~t~ ^i,kn ' ^kn,j (43) 
'^ij.kn — ^i,knikn,j3 ~^ ^iS^kn^ kn,j ~^ ^{^i,kn ^ ^fcn.j) (44) 

_ il-¥^k'S^n'){e-L + ~e^L'r , u)-^ (4^) 

gkn.k'n' - ~ — -r , tt s ^ \^Pk^^Pk'^''-^ d) (45) 

^Vkn ^ Pk'n' ^P^kk' 

with ep^^ = £3 — Epj.^ + 5(Ud — Up). In the following subsections, we will try to estimate these exchange couplings in 
both the strong and weak spin-orbit interaction cases. 
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2. Direct exchange 

Here we require only second order perturbation theory in the direct matrix elements. One obtains the results^^: 



J = 






Ud ' 






4* /^d Id 


Id r-id \ 

ij ji)^ 


■ — 

-L — 




^f^d ~f^d 1 ^{ r^d £~td \\ 



(46) 
(47) 
(48) 



F. Strong spin-orbit interaction 



As discussed above in Sec. IIIID 11 in the strong spin- 
orbit hmit, A ^ |ei,2 — eal, one can obtain effective total 
angular momentum eigenstates with j = 1/2. Choosing 
Eq. ([28]) . and rewriting the corresponding eigenstates in 
the canonical t2n basis, Eq. ([5^ becomes 





2px 


2py 


2p, 


5px 


5py 


5pz 


A, xz 


t 

















A, yz 





t 











~t 


A, xy 











-t 








B, xz 











-t 








B, yz 








t 





-t 





B, xy 


t 


















H31 



1 

1 

71 



) (49) TABLE II: The transfer integrals between the t2g orbitals on 
A and B Ir*^ and the Px,y,z orbitals on the intermediate 
((i)d- + (i — d- ) (50) ions. "2^2," represents the p^ orbital on the 2nd ion in 

t/xz^ i,yz^ i.xyi) , Fig. g] "A, iz" represents the a;2: orbital on the A ion, And the 

entry t on the row of "A, xz" and the column of "2^2," denotes 
the hopping amplitude (transfer integral) from xz orbital at 
A ion to px orbital on 2nd O^^ ion. Other notation can be 
understood likewise. 



in which, we have expressed 0^3^ /a^i in terms of the t2g 
annihilation operator to avoid the position dependence 
of the coefficients. 



1. Superexchange through oxygen ions 

The complicated expression of Eq. (|42|) requires sim- 
plification if we want to have a quantative understanding 
of the exchange coupling. However, some information 
can be immediately obtained from Eq. (j50[) . in partic- 
ular that all ti3,A;„ = 0, which makes J, D^j and F y 
only remain terms with Ci^kn- To simplify further, we 
need some explicit form for the transfer integrals tjm.kn- 
Hence, we will make further approximation that the sur- 
rounding octahedra of Ir*+ are perfect so that we can 
apply the cubic symmetry to find out the nonvanishing 
transfer integrals and also the relation between them, 
which is listed in Table. [TTl for Ir^+ A and B in Fig. g] 
Deviations from these forms should presumably be small, 
since the non-cubic distortion is. 

Based on the transfer integrals listed in Table. |TT1 we 
evaluate the exchange coupling constant J and F ab- 
For bond AB, collecting non-zero coupling constants (ac- 
tually J — 0, Dab =0), we obtain 

Hab = ~JS^aS% + JS\Sl + JS\S% (51) 

with 

J = o I^l'*(2ff2p,,5p^ - 92p^.2p^ - g5p^,5pj ■ (52) 



Since from Eq. (gSl) g2p^,5p^ > 52p,,2p, , g5p^,5p^ , then J > 
0. Thus we find ferromagnetic interaction between the x 
components and antiferromagnetic interactions between 
th y and z components along this link. This corresponds 
to the form in Eq. ^ of the Introduction, with e^j = 
efj = —efj = 1 for this link. 

Because all links are equivalent by point group oper- 
ations, we can deduce the exchange interactions of all 
other bonds by symmetry. The sites A and B correspond 
to point 4 and 8 in our notation in Fig. [31 The result 
is that the exchange interactions on each bond are ferro- 
magnetic between one component, and antiferromagnetic 
between the other two. These principle components are 
always along x, y, or z. We will call a bond in which 
the X component is ferromagnetic a "type x bond" , and 
similarly for y, z. The type of each bond is listed in 
Table. [TTTl This Hamiltonian breaks spin-rotational sym- 
metry strongly. A simple rule can be used to characterize 
the Hamiltonian of a given bond: if bond (ij) is located 
in y-z plane, then the bond is type x bond and has type x 
exchange Hamiltonian; if it is located in x-z plane, then 
the bond is type y bond and has type y exchange Hamil- 
tonian; if it is located in x-y plane, the bond is type z 
bond and has type z exchange Hamiltonian. As a result, 
the three bonds in every triangle (See Fig. ^ have dif- 
ferent exchange Hamiltonian. The ground states of this 
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type X 


type y 


type z 


(1,2) 


(1,3) 


(2,3) 


(3,5) 


(3,4) 


(4,5) 


(5,7) 


(5,6) 


(6,7) 


(4,8) 


(8,9) 


(4,9) 


(8,11) 


(7,11) 


(7,8) 


(1,6) 


(6,12) 


(1, 12) 


(9, 10) 


(2,9) 


(2, 10) 


(10,12) 


(10,11) 


(11,12) 



TABLE III: The bond types of 24 bonds in one unit cell. 
Points and bonds are based on the notation in Fig. [S] "i" is 
used for the points which are simply a translation by a basis 
vector from point "i". 



Hamiltonian will be studied in Sec. UY 



2. Direct exchange 

We consider two Ir atoms A and B, connected by a 
line along the (0,1,-1) direction. There are two prin- 
ciple overlaps. The largest, whose magnitude we denote 
tf, is between the yz orbitals at each atom - this is a 
CT-bond. A secondary overlap, of magnitude ^2, occurs 
between orbitals of the form xy — xz at each site, which 
corresponds to tt bonding. All other overlaps are ex- 
pected to be negligible or zero. This leads remarkably 
to 



Cn' = 0. 



(53) 



The result appears isotropic, despite the strong spin-orbit 
interactions! As a consequence, one obtains only Heisen- 



berg exchange, and 



0! It is remarkable that 



one finds apparent isotropy even though the spin itself 
contains a substantial orbital component. As seen from 
the superexchange calculation above, this is by no means 
guaranteed. 

The first corrections to the strong spin orbit limit are 
linear in the non-cubic splittings, and produce corrections 
to the Heisenberg model. This occurs by a contribution 
to Cjj' of 0{\e3 — ei,2|/A). The leading spin-orbit correc- 
tions to the exchange Hamiltonian are then of the DM 
form, and constrained by symmetry according to consid- 
erations of Sec. nil Al 



G. Weak spin-orbit interaction 

In this part, we are going to look at the weak spin- 
orbit interaction limit, A <C |ei — £2!, |e2 — esl- This is 
the regime which was often studied in literaturep-'^-!^^ 
Standard perturbation treatment can be applied, which 
yields 

-djm'a' (54) 



with timm' introduced previously in Eq. (|24p . Using this 



in Eq. ([33]) reproduces Eq. (|3T|) . 

Keeping the exchange coupling constant to the lin- 
ear order of t — - — r, we can ignore F i,-, as it is of 

2—^31 ' ^ ■" 

0{{ compared with J, thus we only need to 

evaluate J and Di, . 



1. Superexchange through oxygen ions 

Since all the bonds and sites are equivalent, we can 
take ij as bond BA in Fig.Ul Denote the unit directional 
vectors for Di, D2 and D3 as ei, 62 and 63. Ignoring the 
small effect of lattice distortion on these vectors and tak- 
ing the corresponding values for an ideal hyper-kagome 
lattice, we will get ei = ^(0,-1,1), 62 = ^(-1,1,1) 

and 63 = —-^(2, 1, 1). Making the same approximation 
as in previous section, we can evaluate the exchange cou- 
pling constants: 



J 



= 1->BA 



2^ 

71' 



BA 



2A , , 



£2 - £3 

£3 



ei - £3 



ei - £3 



)(55) 



n r> « A I |4^ 3g2px,2px 

V6 £2 - £3 

The three DM components we obtained in Eq. ([55|) are 
not independent from each other. That's because we 
render some symmetry to the system by the transfer 
integrals. Hence, we will still consider all three com- 
ponents to be independent. As discussed in Sec. IIIIB| 
£3 > £2 > £1, additionally, we have gkn,k'n' > 0, then 
we can confer from Eq. ((55|) that J > 0, D2 < 0, is 
probably positive due to a factor of 3 in front of (?2px,2px 
and the smaller denominator of the positive term than 
the negative term, and \Di\ is probably small compared 
to 1 1)2 1 due to the cancellation of positive and negative 
terms. 

Using Eqs. (|55p and ignoring its specific expression, we 
may estimate the strength of DM interactions crudely. 
Then, we estimate crudely 



|A|/J~A/|£i,2-£3|. 



(56) 



Since we assume A <C |£i,2 — £3!; otherwise the perturba- 
tive treatment doesn't holds. We estimate the spin-orbit 
coupling A « 0.4eV, taken from Refs'2Cll)2ll (although the 
reference is not directly relevant to Na4lr30g, we can use 
their spin-orbit coupling as an approximation) . The split- 
ting of the t2g states due to the non-cubic environment 
which determines £1,2— £3 is difficult to estimate. As men- 
tioned in RefsjH, the eg-t2g splitting for [Ir(NH3)( 



IS 

about 5eV. However, if we seek a lower bound on \Di\ we 
can make due with what is probably an over-estimate of 
this splitting. Taking |£i^2 — £3! ~ 5 — lOeV is surely in 
that category, and we therefore find \Di\/J > 0.04 — 0.1. 
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2. Direct exchange 

One can similarly evaluate the induced DM terms at 
first order in the spin-orbit coupling in the case of direct 
exchange. One again obtains a D-vector consistent with 
the symmetry considerations in Sec. 1111 Al 



IV. CLASSICAL GROUND STATES OF THE 
STRONG EXCHANGE ANISOTROPY 
HAMILTONIAN 



In this Section, we will consider the ground states of 
the strongly anisotropic Hamiltonian, Eq. obtained 
in the strong spin-orbit limit from the Ir-O-Ir superex- 
change mechanism. 

Take the triangle A123 in Fig. [3] for example. Bond 
(1,2) is of bond type x; bond (1,3) is of bond type y; 
bond (2,3) is of bond type z. Then for bond (2,3), the 
Hamiltonian is 



n 



(2:3) 
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Clearly 7^2,3 is minimized if 



TABLE IV: The basis vectors for the ground state spin con- 
figurations of the strong spin-orbit Hamiltonian. 



CLASSICAL GROUND STATES INDUCED 
BY DZYALOSHINSKII-MORIYA 
INTERACTIONS 

A. Order due to second component only 



02 02 

a^,y _ _a^,y 

O2 — ^3 



(58) 
(59) 



In general, for each bond, the energy is minimized if 
the ferromagnetically interacting components of the two 
spins involved are parallel, and the antiferromagnetically 
interacting components are antiparallel. We can search 
for unfrustrated ground states by demanding this on ev- 
ery bond. Fixing one spin, its neighbors are therefore 
determined, and from them further neighbors, etc. It is 
straightforward to verify that in this procedure no con- 
tradictions are encountered despite the presence of loops 
on the lattice. In this way all classical ground states are 
determined from the choice of a single initial spin. Thus 
the Hamiltonian is unfrustrated, and we have found its 
full set of classical ground states. Mathematically, we 
can write the full spin configuration as 



s, = s^Vi. + syv2. + s'V3 



3,i 



(60) 



where s = (s^, s^, s^) is a unit vector, and Va^i is the vec- 
tor Va corresponding to the ith spin in Table. ITVl We see 
that the ground states are parametrized by two continu- 
ous parameters - the angles specifying the orientation of 
the initial spin, or of s. This is actually an accidental de- 
generacy, since the system has only discrete space-group 
symmetries, but it is very small. Still, it should be re- 
duced to a discrete degeneracy by perturbations such as 
quantum or thermal fluctuations, or additional interac- 
tions, which will select a subset of these states. 



In Sec. mil we found that the direction of the DM vec- 
tor for a single bond is arbitrary, i.e. not determined 
from symmetry considerations, and not calculable from 
microscopic theory without a more detailed understand- 
ing of matrix elements than we have at present. A gen- 
eral solution for the ground state with such an arbitrary 
DM vector is quite difficult, because different triangles in 
the hyper-kagome lattice are located in different planes. 
In this subsection, we will consider the special case in 
which the DM vector is normal to the triangular plane, 
i.e. Di = D3 = (see Fig. [3]). This is a helpful start- 
ing point for the more general case which we will address 
thereafter. 

As in the case of the nearest-neighbor kagome antifer- 
romagnet, a nonvanishing D2 (here, by D2 we mean the 
component of DM vector which is normal to the kagome 
plane) selects coplanar ground states with 120° spin ori- 
entations on each trianglej^^ii^ii^ These are the only con- 
figurations in which the Heisenberg interactions on a tri- 
angle are minimized (i.e. the total sum of spins on a 
triangle is zero) and the DM interaction is minimized at 
the same time. In the kagome lattice, however, the copla- 
nar ground state manifold is highly degenerate, since ro- 
tating the spins in a single hexagon about the normal 
axis of the kagome plane by arbitrary angle generates a 
new ground state from any other one. In contrast, for 
the hyper-kagome lattice, the non-coplanar nature of dif- 
ferent triangles reduces the degeneracy to just a pair of 
Kramer's degenerate (reversed) states. One of them is 
drawn in Fig. [21 the other one is generated by reversing 
all the spin directions. The chirality of the hyper-kagome 
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A 




FIG. 6: (Color online) Illustration of spin direction at shared 
corner of two neighboring triangles. Black arrows indicate the 
DM interaction path. B, C and E are in the same line. 



lattice makes this state a ground state only for D2 < 0. 
With the other sign of D2 , the DM and Heisenberg inter- 
actions cannot be simuhaneously satisfied. We wiU call 
these states uncanted "windmill" states - see Fig. [2] 

To see that the uncanted windmill states are the only 
classical ground states, see Fig. [6l Starting from triangle 
ABS - denoted AABC - a nonvanishing D2 component 
prefers a coplanar spin configuration, which requires that 
spin A, B, C should lie in the AABC plane and at 120 an- 
gles as dictated by the Heisenberg interaction. The same 
applies to ACDE. However, AABC and ACDE are not 
in the same plane, which confines the spin orientation of 
site C to be aligned with the intersection line of AABC 
and ACDE. We apply this result to all spins, and the 
Heisenberg interaction will select two states, which si- 
multaneously minimize the DM interaction with D2 < 0. 
The magnetic unit cell of the windmill state is the same 
as the chemical cell. 

The result that D2 < completely removes the mas- 
sive but accidental ground state degeneracy of the hyper- 
kagome is quite dramatic. A classical antiferromagnet 
with this interaction will clearly order at low temper- 
ature, and the drastic reduction in degeneracy suggests 
that even for a quantum system, the suppression of quan- 
tum fluctuations by D2 may be large. Before turning to 
this, we continue with the analysis of classical ordering 
in the remainder of this section. 



B. Magnetic representational analysis of space 
group 

Representational analysis of the magnetic space group 
has proven to be a useful tool to extract informa- 
tion about low temperature ordered phases using lattice 
symmetryj^i2ii2^i^5 . The idea is to consider those types 
of magnetic order which can be reached by a continuous 
transition from the paramagnetic state, which has the 
full space group symmetry. Though there is no a pri- 
ori reason why the ground state configuration need be of 
this type, this is a convenient way to generate candidate 
magnetically ordered states. In principle, one may iter- 
ate this procedure to generate lower temperature ordered 
states, generating all possible ordered phases. 
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TABLE V: The basis vectors of one dimensional irreducible 
group representations of the space group P4i32 appearing in 
the magnetic representation with k — (0, 0, 0). 



The operators of the space group act on both the po- 
sition of the magnetic ion and on the components of the 
spin vectors. The combination of these two results are de- 
scribed by the magnetic representation T. The magnetic 
representation for a particular site can be decomposed 
into contributions from the irreducible representations of 
the little group 

r = 5]n^r^. (61) 

For Na4lr308, the space group is P4i32 (although it 
can also be P4332, the results should be equivalent),— 
and the Ir*+ ions sit on the 12c?. Here, we only focus the 
simplest case when the propogation vector k = (0,0,0). 
A program called "SARAh"— is used to do the decom- 
position of magnetic representation 

r = ir^^' + 2r^^' + sr^^' + srf + srf ^ , (62) 

in which, the superindex represents the dimension of the 
irreducible representations, and the subindex counts the 
irreducible representation. 

Landau theory requires that only one representation 
can be involved in a critical transition, and so with this 
constraint there are only five possible magnetic struc- 
ture for k = (0,0,0). Even within this decomposition 
and Landau theory constraints, for certain representa- 
tions (r3,r4,r5), there still remain a lot of degrees of 
freedom because of the multiple basis elements in these 
3-dimensional representations. For simplicity, we only 
discuss the one dimensional representation Fi, The 
basis vectors for these two representations calculated are 
given in Table. |Vl 




FIG. 7: (Color online) The basis vector tps in Table. El 



The physical interpretation of these representations is 
as follows. The basis vector ipi is nothing but the C2 
rotation axis at every magnetic ion. The basis vector 'ip2 
gives the spin directions of classical uncanted windmill 
state discussed above (see Fig. [5]). The third basis vec- 
tor -03 may be obtained as the axis which is normal to 
both C2 axis and the spin direction in ^^2 (see Fig. [2|) . 
Note that these three basis vectors at each site form an 
orthonormal basis for the spin coordinates. 

Evidently r2 is related to the DM interaction, at least 
to the D2 component. But what about Di and D^l Let 
us consider following situation. Starting from an ordered 
ground state with D2 < 0, we turn on an infinitesimal 
Di or D3 component. The spin at site i can be written 
as 



S, = Jl-(ai)2-(a3)2e2 + aiei + a?e: 



3-3 
i 1 



(63) 



where ei, 62 and 63 are simply the three orthogonal unit 
vectors given by basis ipi^ ^"2 and ip3, and a} and af are 
small corrections to the ordered ground state due to the 
introduction of an infinitesimal Di (or D^) component. 
We plug Eq. [521 into the Hamiltonian, and expand to the 
2nd order in a], af and Di (or Z?3).To linear order 
vanishes. Thus the ground state spin configuration with 
negative D2 and an infinitesimal Di or D3 component 
is related to 'ip2 and -03 • The irreducible representation 
T^2^ is relevant to the magnetic structure when the DM 
interaction is present. 

Now, we proceed by assuming that the ground state 
configuration ip of the more general case, when Di, D2 
and arc all present in the system, is a linear superpo- 
sition of basis vectors ip2 and ^3 : 



Si = cosxt/jj i + sinxip 



3,i 



(64) 



where xp^ ^ is the vector ipa corresponding to the i^^ spin 
in Table |Vl Evaluating Eq. ^ for spin configurations of 



FIG. 8; (Color online)The spin configuration of site 
1, 2, 3, 4, 5. The dashed blue lines are the canting axis of cor- 
responding spin. 



this form gives 

n/N = 2{-3V2Di + 5VSD2 - V6D3 - 3J) 

+2^/3[(\/6L>l + D2 + V2D3 - VSJ) cos {2x) 
-{VSDi - 3D3) sin (2a;)] , (65) 

where N is the number of unit cells in the lattice, not the 
number of spins (which is equal to 12iV). Minimizing the 
Hamiltonian with resepct to a;, we can find the canting 
angle x is given by 



cos (2a;) 
sin (2a;) 



V3J + V6D1 +D2 + V2D3 



W 



W 



(66) 
(67) 



where we have defined 



W = y{V6Di + D2 + V2D3 - V3J)2 + 3{Di - ViDzf 

(68) 

for convenience. 

Fig. [51 is an example of this canted state when Di — 
O.IJ, D2 = -0.04 J and D3 = 0. We only plotted the spin 
configuration of A123 and A345 in Fig. [31 Because these 
states are obtained by smoothly introducing a ^3 compo- 
nent into the uncanted windmill states found in Sec. |Vl 
we will call this state a canted windmill state. The canted 
moment disappears not only when Di = D3 — Q, but also 
for Di — ^/SDs, at which point it degenerates into the 
uncanted windmill state. Regardless, it is also interest- 
ing to note that because Di, D3 <C J, from Eq. ((67l) the 
canted moment is much smaller than the coplanar com- 
ponent, which indicates the dominance of the coplanar 
spin configurations. Similar features have been found in 
other studies. 

In the above treatment, we have assumed the ansatz in 
Eq. (|64p , which is not guaranteed to find the global mini- 
mum energy state, and gives no guidance as to where this 
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form of the ground state breaks down. In next section, 
we will consider this question from a different point of 
view. 



C. Mean field spherical model 

In this subsection we approach the general problem of 
finding classical ground states of the Hamiltonian from a 
different point of view. The problem is difficult because 
in addition to minimizing 7i, which is quadratic in spins, 
we must also satisfy constraints that each spin have fixed 
magnitude |Si| — 1. The large number (equal to the 
number of spins) of these constraints makes what other- 
wise would be a simple quadratic minimization problem 
difhcult. Here we replace these many constraints by a 
single one, 

^|S,f = 127V, (69) 

i 

where, as elsewhere in the text, we define N as the num- 
ber of unit cells for convenience. This is the "spherical 
model" , and is exactly soluble at both zero and non- 
zero temperature. At zero temperature, the spherical ap- 
proximation must give a lower bound to the true ground 
state energy, since minimization is conducted with less 
constraints than in the physical spin model. Because of 
this observation, this approach can indeed often be used 
to construct physical ground states. This "Luttinger- 
Tisza" methodi ^^'^^1^^ consists of finding a subset of 
ground states of the spherical model which respect the 
spin normalization constraints of the physical problem. 
Any such states must be ground states of the full Hamil- 
tonian. Moreover, when such states exist, they exhaust 
the full set of physical ground states. However, it is not 
always possible to find any ground states of the spherical 
model which satisfy the normalization constraint. If not, 
it simply means that ground state energy of the physi- 
cal problem is strictly larger than that of the spherical 
model, and the Luttinger-Tisza method fails. Generally, 
the Luttinger-Tisza construction is less effective on lat- 
tices with a large number of sites in their basis. For the 
hyper-kagome lattice with a 12 site basis, our expecta- 
tions should not too high! Nevertheless, in some range of 
phase space, we will indeed find physical ground states 
from this approach. More generally, at non-zero temper- 
ature, the spherical model may be a useful approximation 
even when it fails to produce exact ground states at zero 
temperature. 

Minimizing the quadratic Hamilton in Eq. ^ with the 
single global constraint in Eq. (j69p is a standard prob- 
lem, which is solved by finding the eigenvectors of the 
Hamiltonian matrix (coefficients of the quadratic form of 
spin components) with minimum eigenvalues. By trans- 
lational invariance, the eigenfunctions have the Bloch 
form, i.e. are quasimomentum eigenstates. Hence it is 



useful to Fourier transform Eq. ([2]): 

n = Nj2J2J2^UW.^(^)Qji^) ■ (70) 

k i,j 

Here 

S^RrO = ^Qr(k)exp(ik.R„) (71) 

k 

with V, fi index of spin vector components, R„ is the po- 
sition of unit cell, i and j are the sublattice index and 
is the Fourier transformed Hamiltonian matrix in 
the Bloch representation (which is 12 x 3 = 36 dimen- 
sional because of the multiple basis sites and spin com- 
ponents). We need to minimize Eq. ([70)1 subject to the 
soft constraint Eq. ([7^ . which can be expressed as 

12 

^^Ql(k).Q,(k) = 12. (72) 

i=l k 

Minimization is equivalent to find the minimum eigen- 
values (and corresponding eigenvectors) of L^^^(k). We 
did this numerically for every k, and found the global 
minimum for every (Di,D2), {D2,D'i) and {Di,D^) 
pairs. With this approach, phase diagrams in D1-D2, 
D2-D3, D1-D3 parameter spaces have been plotted in 
Fig. [H In a wide region of the phase diagrams, the min- 
imum eigenvalue is realized for k = (0,0,0). In this 
case, the corresponding eigenfunction can be chosen to 
satisfy the normalization constraint on every site, and 
so an exact ground state is found. This ordered state 
is in fact precisely the canted/uncanted windmill states 
we proposed in previous section. Thus in the regions 
for which k = (0,0,0) is indicated in the figures, this 
analysis proves that these windmill states are the exact 
global ground states. In a considerable large regions of 
the parameter spaces, we get canted/uncanted "wind- 
mill" states. 

In other broad regions of the phase diagram, the spher- 
ical model predicts ordered states with incommensurate 
wavevectors, i.e. in which k has irrational projection 
onto reciprocal lattice vectors. This is indicated simply 
as "incommensurate phase" in the figures. In most cases 
we have studied, the incommensurate wavevectors are 
located around 0.85(7r, tt, tt) and its eight equivalent mo- 
menta 0.85(±7r, ±7r, ±7r). However, in this region of the 
phase diagram, we are unable to construct a linear com- 
bination of eigenfunctions which satisfies the local con- 
straint on the spin magnitudes. Thus the incommensu- 
rate ground state of the spherical model does not immedi- 
ately imply a corresponding ground state of the physical 
model. It is possible that the region of phase space occu- 
pied by the windmill states is actually expanded beyond 
what is shown here by this effect. Most likely, ground 
states with large unit cells or incommensurate order do 
exist in the physical model, but are more complex than 
those of the spherical approximation, and with some- 
what higher energy. Even in the spherical model, we 
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FIG. 9: (Color online) phase diagram in D1-D2, D2-D3 and 
D1-D3 parameter spaces. The uninvolved DM vector com- 
ponent in each figure is set to be 0. The red reference lines 
(axes) are not phase boundary. 



see that in the incommensurate region, while the DM 
interaction removes much of the frustration-induced de- 
generacy, the enlargement of the unit cell implies a larger 
residual ground state degeneracy, and hence less effective 
removal of frustration than in the k = (0, 0, 0) regions. 



VI. QUANTUM EFFECTS 

A. Numerically constructed Bogoliubov 
transformation 

In previous sections, the spins were treated classically 
and classical ground states were obtained. In this sec- 
tion, we discuss the quantum effect in the formalism of 
linear spin wave theory. In certain regions of D1-D2, 
D2-D^ and D1-D3 parameter space, we have an ordered 
ground state. We will use the Holstein-Primakoff Boson 
approach to explore the quantum effectsi^ Now express 
the spin operator as follows. 



Si(Rn) • Si 



ord 



Si(R„) • 



V2S 



Si(Rn) • {S° 



el) 
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J2S 



a»(R„) (73) 
(a,(R„)-a|(R„))(74) 

(a,(R„) + a|(R„)) (75) 



where S°^'^ is the unit vector along the spin order, b\ is 
the C2 rotational axis at site i introduced in Sec. IV Bl 
a\ and are the creation and annihilation operators of 
Holstein-Primakoff bosons at ith sublattice of unit cell 
at position R„, and we only keep the lowest order of a\ 
and ai. Under this transformation, the Hamiltonian can 
be written as 

S 



n 



H 



- ^ A,,(k)al(k)a,(k) -I- B,,{k)al{k)a]{-k) 

(76) 



where we have dropped the constant term and high order 
terms. Here, Aij(k) and Bij^k) are the coefficient ma- 
trix we end up with after doing Fourier's transform on 
the creation and annihilation operators. The Fourier's 
transform we used is 

ai(R„) ^ ^ Oi exp (-ik • R„) . (78) 

Since there are 12 sublattices, using the analytical Bo- 
goliubov transformation is hopeless to diagonalize the 
Hamiltonian. Here, we will use a numerically constructed 
Bogoliubov transformati on (NC BT) introduced and dis- 
cussed in detail by Ref. I3lll32l to diagonahze Eq. (175]) . 
find the spin wave energy gap and calculate the quantum 
corrections to the classical order. Write Eq. (|76p as 



7i = ^Xt(k)H(k)X(k) 



(79) 



where 
X(k) 

H(k) 



(ai(k) . . . ai2(k), a\ (-k) . . . a\,{-k)f (80) 

S( A{k) B{k) \ 
2 \B*{~k) A*i~k) ' 



(81) 
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and the hermiticity of 7i requires that 

A,,(k) = A^(k) (82) 
By(k) = B,,(k) . (83) 

We now introduce the canonical transformation 

X(k) = Q(k)Y(k) , (84) 

where, Y(k) is given by 

Y(k) = (foi(k) . . . 6i2(k), b\{-k) . . . b\,{-k))^ , (85) 

and satisfies 

[6,(k),6](k')] =<5,,5k^k' . (86) 

The transformation Q is required to to diagonahze the 
Hamihonian as 

Qt(k)H(k)Q(k) =A(k) , (87) 

where A(k) is the diagonal eigenvalue matrix 
whose diagonal matrix elements are given by 
(ei(k),- • • ,ei2(k),ei(-k),- • • ,ei2(-k)). Using this 
transformation, the quantum correction to the classical 
spin polarization can be written as 

= I^E(«l(k)«»(k)) 

= ^y"(Xtx) - - , (88) 

k 

At zero temperature, further making use of Eq. ([M]) . 
Eg. ([55)1 can be expressed as 

k i=l 

If we find the canonical transformation Q(k), the energy 
spectrum can also be obtained. With the energy spec- 
trum, we can find the spin wave energy gap, A. Some 
care must be taken as the numerical construction of the 
Bogoliubov transformation is effective only when there is 
an energy gap. 

B. Quantum corrections and spin wave gaps 

By the method described in last section, we carried out 
the numerical procedure described in Ref. 31,32. Taking 
spin iS* = i, we numerically construct the Bogoliubov 
transformation for every k, and find its contribution to 
zero temperature quantum correction, dS^ and energy 
levels at every k to extract the spin wave gaps. The nu- 
merical results are plotted in Fig. (TUl Fig. [Til Fig. [T^] 
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FIG. 10: (Color online) The dependence of quantum correc- 
tions and spin-wave gaps on the DM vector components. In 
the two figures, we set Di = Ds = and vary D2 . 60 x 60 x 60 
momentum points have been used to generated the data. 
No change has been found in quantum corrections and gaps 
within computer resolution compared with 50 x 50 x 50 mo- 
mentum points. 



and Fig. [131 Corrections dS larger than 1/2 have been 
truncated to 1/2. In these figures, two components of 
the DM vector are kept constant while the third is var- 
ied. The ordered regions (the third varying DM vector 
component) of these figures can be found the Fig. [9l 

In these figures, spin wave gap is non- vanishing, so our 
calculation is valid. It is easy to find the quantum be- 
havior also resembles the classical one: the different DM 
vector components have different effects in quantum cor- 
rections, which is similar to the effect of DM vector com- 
ponents in favoring canted "windmill" state in Sec. IV Bl 
In the DM magnitude studied in these figures, the quan- 
tum corrections are pretty large. Even in the case when 
D2 — — 0.09J and Di ~ D3 — 0, the quantum correction 
is about 50%. 

As a general rule, one observes that the quantum cor- 
rections decrease steadily as one goes deeper into the 
k = classically ordered region. If we crudely suppose 
that dS > 1/2 is indicative of the destruction of order 
by quantum fluctuations, we may expect broad regions 
of quantum spin liquid states occurring in and near the 
incommensurate regions of the classical phase diagram. 
This range of DM vectors then may be possible candi- 
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FIG. 11: (Color online) The dependence of quantum correc- 
tions and spin-wave gaps on the DM vector components. In 
the two figures, we set D3 = and vary Di with two fixed 
D2 values (D2 = -0.08J(in blue) and D2 = -Q.05J(in red)). 
16 X 16 X 16 momentum points have been used to generated 
the data. No change has been found in quantum corrections 
and gaps within 1% compared with 10 x 10 x 10 momentum 
points (Same for Fig. [12] and Fig. [T3j. 

dates for application to Na4lr308- 

C. Comparison with exact diagonalization 

In order to partially confirm our results in last sec- 
tion, we performed numerical exact diagonalization for 
5 = 1/2 spinsi^^i^ We took six triangles with thirteen 
sites and used a Heisenberg model plus DM interaction 
with only D2 ^ 0- We plot the resulting specific heat in 
Fig. [m The gap in each case can be inferred from the 
plot by the temperature below which the specific heat 
becomes negligible. As we found in previous section, the 
more negative D2 is, the greater the gap will be. At 
low temperatures in Fig. 1141 the sequence of the curves 
agrees with what they should behave according to spin 
wave gaps in Fig. [TOl 

Similarly, we also look at the case when only D3 com- 
ponent is present by taking = ±0.10 J, ±0.25 J. Ac- 
cording to Fig. [131 the spin wave gaps of D3 = ±0.10 J 
are close to each other, and the spin wave gaps of 
D3 = ±0.25J are also close to each other, but much 
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FIG. 12: (Color online) The dependence of quantum correc- 
tions and spin-wave gaps on the DM vector components. In 
the two figures, we set D2 = 0, D3 = —0.1 and vary Di. We 
also did the same thing with D2 = 0, D3 = 0.1 and varying 
Di, and the quantum corrections always break the classical 
order completely. 

larger than the previous cases. In Fig. [13 we see that 
both curves of D3 = ±0.10J and D3 = ±0.25J nearly 
overlap at low temperatures, and their sequence agrees 
with the magnitudes of the spin- wave gaps. 



VII. DISCUSSION 

In this paper we have studied the effect of spin-orbit 
interactions in the hyper-kagome lattice of Na4lr308. A 
crucial physical parameter is the strength of atomic spin- 
orbit coupling relative to non-cubic crystal field split- 
tings. In the strong spin-orbit limit, Ir-O-Ir superex- 
change produces highly anisotropic effective spin inter- 
actions, with 2/3 antiferromagnetic and 1/3 ferromag- 
netic couplings between spin components. This Hamil- 
tonian turns out to be unfrustrated, and has a small set 
of classical ground states. We speculated this even the 
S — 1/2 quantum model is likely ordered with a substan- 
tial critical temperature, inconsistent with experiments 
on Na4lr308. 

By contrast, strong spin-orbit interactions and direct 
Ir-Ir exchange, or weak spin orbit interactions, all induce 
an Heisenberg-like Hamiltonian with a small correction 
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FIG. 13: (Color online) The dependence of quantum correc- 
tions and spin-wave gaps on the DM vector components. In 
the two figures, we set Di = and and vary D3 with two D2 
values {D2 = 0(in blue) and D2 = -0.04J(in red)). We also 
did the same thing with D2 = 0, -D3 = 0.1 and varying D\, 
and the quantum corrections always break the classical order 
completely. 
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FIG. 14: (Color onhne) The specific heat of six triangle with 
negative D2 component DM interaction. Along the thin ver- 
tical line, from top to bottom D2 value of each curve increases 
from -0.09J to -O.OIJ with a step 0.01 J. 
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FIG. 15: (Color online) The specific heat of six triangle with 
Dz component DM interaction. Along the thin vertical line, 
the upper curves have ID3I = 0.25 J, the down curves have 
ID3I =0.10J. 



of the form of a Dzyaloshinskii-Moriya (DM) term. The 
orientation of the DM vector, which is not determined 
by symmetry or our microscopic considerations, deter- 
mines the extent to which the frustration of the Heisen- 
berg model is relieved. In one region of phase space, 
frustration is fully relieved, and the DM interaction stabi- 
lizes a magnetically ordered "windmill" state, with nearly 
but generically not quite coplanar moments. Quantum 
fluctuations, which we assessed by spin wave theory, are 
sufficiently suppressed that we may expect this order to 
persist even for spin-1/2 spins (as in Na4lr308) in part of 
this region. In the remainder of phase space, the frustra- 
tion is not fully removed, and the classical ground states 
break the lattice periodicity and may be incommensu- 
rate. We argued that in this regime, the classical ordering 
is fragile and may be destroyed by quantum fluctuations 
for S* = 1/2 spins. 



A. Zero temperature susceptibility in quantum 
spin liquids 

Part of the motivation of the present study was the ob- 
servation in Na4lr308 that the susceptibility x tends to a 
constant at low temperature, despite the approximately 
quadratic decrease of specific heat. We argued that this 
combination, which implies a diverging Wilson ratio as 
T ^ 0, is likely indication of spin-orbit interactions. In- 
deed, on general grounds, a constant zero temperature 
susceptibility is expected when SU(2) spin-rotation sym- 
metry is broken. The situation of weak DM interaction 
is quite common in frustrated magnets, and may allow 
this behavior quite broadly. Therefore it is interesting to 
consider more generally how this occurs in the presence 
of weak DM coupling. We have not so far addressed the 
magnitude of this zero temperature susceptibility. 

Presuming the DM interaction to be relatively weak, 
the magnitude of x should be understood in terms of the 
correlations the spins would have in the underlying sys- 
tem without DM. Various SU(2) invariant phases lead to 
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rather different behaviors. Generally speaking, one ex- 
pects the most suppressed x for systems with the least 
low-energy spin fluctuations in the absence of DM. Prob- 
ably the most extreme example is a Valence Bond Solid 
(VBS) or dimer state, in which the eigenstatcs can be 
approximated by those of a single partition of the sites 
into pairs of spins which are coupled to each other only 
within the pairs. Such a VBS phase has a gap of order 
J to all excitations, including the elementary triplets. A 
simple calculation by second order perturbation theory 
of the susceptibility shows that it is indeed non-zero, and 
of order 
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One may also estimate the magnitude of x for various 
phenomenological gapless spin liquid ground states per- 
turbed by DM. The general arguments follow scaling 
theory. We presume the gapless spin liquid is a criti- 
cal phase in the renormalization group sense, described 
by a scale invariant field theory. Introduction of DM 
interactions breaks SU(2) symmetry, and allows oper- 
ators Oa breaking SU(2) to be added to the effective 
action/Hamiltonian. Generically, these appear with co- 
efficients proportional to D. In the simplest situation, 
there is a single such operator Oa with the smallest scal- 
ing dimension A. In most cases of interest, we expect 
A < d + z, where d is the spatial dimension and z is 
the dynamical critical exponent (z = 1 is common). In 
this case, the presence of this operator in the Hamilto- 
nian constitutes a relevant perturbation. Then, if the 
susceptibility a,t D — Q behaves as x ~ T^, we expect 
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where /i, v are spin components x,y, z. The operator Oa 
is expected to break SU(2) down to some subgroup. This 
may contain either one or zero residual U(l) spin rotation 
axes. The susceptibility normal to this axis, if it exists, 
is expected to be constant at low temperature. If no such 
axis exists, then the susceptibility will be constant in all 
directions. In cither case, we must have 



for \X\ > 1. 
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Here A^,^ is a symmetric tensor with either 2 or 3 non- 
zero eigenvalues, in the cases with one or zero residual 
U(l) symmetries, respectively. One thereby obtains 



X^.AT = 0)-\D\^T7^A^,. 
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As an example, consider the 2d "Dirac" spin liquid 
with point nodes on the kagome lattice studied by Her- 
mele et aJ^. There, the dominant operator indeed pre- 
serves a single residual U(l) symmetry. Its scaling di- 
mension IS estimated as A w 2 - 32/(87r2) « 1.6 (based 
on a calculation for a generalized model with large num- 
ber, Nf, of flavors of fermions, evaluated for the physical 



case Nf = 4). Taking d = 2,z = 1, P = 1 as appropriate 
for this case, we find, restoring units 
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Here XJ_ is the susceptibility in the x-y plane perpendic- 
ular to the conserved U(l) spin axis. We see that the 
dependence on D is sub-linear, making for a very large 
susceptibility even for rather small D/ J. 

It is noteworthy that the scaling prediction above 
should obtain regardless of the other properties of the 
system in the presence of DM interaction. The relevance 
of Oa at the spin liquid fixed point indeed implies that it 
drives the system into a different phase, which may not 
be a spin liquid at all. This is believed to be the case for 
the above Dirac spin liquid, for which the resulting state 
is expected to be magnetically orderedi^ 



B. Other possibilities 

One may wonder whether the weak and strong spin 
orbit limits are the only possibilities for Na4lr308, and 
whether they might be distinguished more directly. Prob- 
ably the principal difference in the two cases is the sign 
of the g-f actor. In the weak spin orbit limit, one has ap- 
proximately M w — 2/iBS, while in the strong case, we 
found M w +2^bS- While these lead to identical Curie 
laws, they are physically distinct (note that one cannot 
reverse the sign of S and maintain its canonical commu- 
tation relations). It should be measurable in other exper- 
iments such as nuclear magnetic resonance. Microscopic 
reasoning gives no reason why the Ir"*"*" spins might not 
be in an intermediate situation between the two extreme 
limits. However, in this case one would expect a (7-factor 
in between these two values, i.e. with substantially re- 
duced magnitude. A large deviation would seem to be 
inconsistent with the measured spin susceptibility. 

This tends to support the notion that Na4lr308 is 
either in the strong or weak spin orbit limit, and not 
in between. Given the incompatibility of the strongly 
anisotropic Ir-O-Ir superexchange Hamiltonian in the 
strong spin-orbit case with experiment, we are led to 
believe the weakly anisotropic Hamiltonian with DM 
interactions is most appropriate (we note that "weak 
anisotropy" still allows for \D\/ J ~ 0.1 which would have 
strong effects on the low energy physics). This, however, 
still leaves open the issue of weak versus strong spin-orbit 
interactions. Though susceptibility experiments do not 
distinguish the two cases, they are physically distinct, 
and could be discriminated by magnetic resonance meth- 
ods, for instance. So far as we are aware, all prior mea- 
surements of ions capable of this distinction have 
been interpreted in terms of the strong spin-orbit sce- 
nario (see for example, ReffSGl). This fundamental phys- 
ical question in Na4lr308 warrants further investigation. 

Could there be another scenario? We cannot rule out 
the possibility that other interactions might play a role. 
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Perhaps further neighbor exchange or spin-lattice cou- 
pling might be significant. These are important subjects 
for future theoretical studies. 
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